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K. C. SHIN 

Abstract. We consider the non-Hermitian Hamiltonian H = — ^ + P(cc^) — (ix)^""*"^ 
on the real hne, where P{x) is a polynomial of degree at most n > 1 with all nonnegative 
real coefficients (possibly P = 0). It is proved that the eigenvalues A must be in the sector 

I ^'^zM — 2n+3 - ^^^^ ^'^^ ^ ~ ~lx^ ~ {^^)^ I establish a zero- free region of 

' the eigenfunction u and its derivative u' and we find some other interesting properties of 

, eigenfunctions. 

O 



in 



Preprint. 
1. Introduction 



> 

^ ■ We are considering the eigenproblem 

O ■ 

■ (1) -u"{x) + [P{x^) - {ixf'"'^^]u{x) = \u{x) for -oo < x < oo 

o ; 

^ ■ with m(±oo) = 0, where P{x) is a polynomial of degree at most n > 1 with all nonnegative 



real coefficients (possibly P = 0). 

This is an example of a class of problems, the so-called PT-symmetric non-Hermitian 
Hamiltonian problems, which have arisen in recent years in a number of physical contexts 



1^ , |T5|, D. Bessis conjectured in 1995 that: 



Conjecture. Eigenvalues of H = — — (ix)^ are all real and positive. 

■ Many numerical and asymptotic results [^, |^, |^, |^ support this conjecture. And later 

for n > 1 it was conjectured that the equation ([l|) also has positive real eigenvalues, under 
different boundary conditions p|. However, there is no rigorous proof of this to date. 

This paper is organized as follows: In Section ^ we prove that eigenvalues of the equation 
(0) lie in the sector | argA| < 2^^- This goes part way to proving that the eigenvalues are 
real and positive. We generalize this result to H = — ^ + [P(x^) + ixQi^x"^)] for some real 
polynomials P and Q. In particular, for the potentials —{ix)^ and -|- igx^ with any real g, 
we have that | argA| < |. Then next in Section ^, for the case H = — ^ — {ixY, we fairly 
precisely locate the zeros of the eigenfunctions and their first derivatives in the complex 
plane. Conversely we find a large zero-free region. In Section ^, still with H = — ^ — (ix)^, 
we find a large class of polynomials that are orthogonal to |mP on each horizontal line. And 
finally in the last section, we discuss related open problems. 
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For the rest of Introduction, we provide some more background information on ([l|). First, 
a PT-symmetric Hamiltonian is a Hamiltonian which is invariant under the product of the 
parity operation V{: x ^-^ —x) and the time reversal operation T(: z h-* —i). Certainly (|I|) is 
PT-symmetric while, for example, —-^+x—{ixY is not PT-symmetric. If if = —-^ + V{x) 
is PT-symmetric, then V{—x) = V{x) and so ReV^(x) is an even function and ImV^(a;) is 
an odd function. Hence if V{x) is a polynomial, then V{x) = P(x^) +ixQ{x'^) for some real 
polynomials P and Q. 

Next by the work of Caliceti et al. [T^, it is known that the PT-symmetric Hamiltonian 



H = — ^ + — g{ixY has discrete spectrum, for g real, and these eigenvalues are positive 
real if g is small enough. However, there are some PT-symmetric Hamiltonians that have 
no eigenvalues |jl8|, §1], or non-real eigenvalues [jl2|, footnote on page 26]. 

Lastly, for any A G C there are two linearly independent solutions of (|l]), if the boundary 
conditions are not imposed. In generic cases, the solutions blow up at both +oo and — oo, 
while in exceptional cases, the solutions decay to zero as x approaches +oo or —oo. Only 
in very exceptional cases (when A is an eigenvalue!) does one find a solution that decays to 
zero at both +oo and — oo (see Lemma |l] for details). 



2. The eigenvalues lie in a sector 

In this section, we prove that the eigenvalues A of (|l|) lie in the sector | argA| < and 
we extend this result for more general cases. To do this we will use results of Hille [0, §7.4]. 

For any A G C the equation (|1|) without the boundary conditions allows two linearly 
independent solutions. If u{x) solves the ODE (0), then since Piz^) — {izY^^^ is an entire 
function (analytic in the whole complex plane), there exists an entire function u{z) which 
agrees with u{x) on the real line and satisfies —u"{z) + [P^z"^) — {izY^^^]u{z) = Xu{z). We 
begin by describing the asymptotic behavior of u near infinity. Recall that degP < n. 



Definition. Let 



j arg(i2"+i) 

Ui = Zn 

^ 2n + 3 2n + 3 



2n+3 
2n+3 



if n is even, 
if n is odd. 



We define Stokes regions 



Sj = {z e C : 9j < aigz < 9j+i}, 



for j = 0, 1, 2, 2n + 2. And for notational convenience, we define Sj+2n+3 
Also we denote 

Sj^^ = {z E C : 9j + e < aigz < 9j^i — e}. 



Sj for all j. 



for < e < 



2n+3 ■ 
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Figure 1. For n = 1; the solid line is the real axis and the dotted rays are 
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the critical rays, argz = Qj = j^, |, , ^ and ^'^'^ 



Notice 6j is neither nor it. Thus the negative and the positive real axes lie within two of 
the Stokes regions (see Figure 1). We call these the left- and the right-hand Stokes regions, 
respectively. Also we call the rays {argz = 9j} "critical rays". 



Lemma 1. Every solution of —u"{z) + [P{z'^) — {izY^^^]u{z) = Xu{z) is asymptotic to 

2n+l 

(2) {const. )z ■* exp 



M — (1 + 0(1)) 



"2n + 3 

as z ^ oo in Sj^^, for each < e < g;^- The error o(l) is uniform in argz in the sense that 
lim^^oo sup{|o(l)| : z G Sj^^, \z\ = r} = 0. 

^4/50 u has infinitely many zeros in C but only finitely many in UjSj^e, for each < e < 



2n+3 ■ 

The asymptotic expressions imply in particular that in each Stokes region, u{z) either 
decays to or blows up, as z approaches infinity in Sj^e- 

Proof. See Hille's book §7.4] for a proof of a more general result. An outline of the proof 
is as follows: Hille first transforms the equation into another complex Z-plane by using 
the Liouville transform. And then he compares u with the solutions of the sine equation 
w"{Z) +w{Z) = and finally transforms back to the original complex 2;-plane. So the above 
asymptotic expressions are the asymptotic expressions for solutions of the sine equation 
(in the Z-variable) expressed in terms of the original z-variable. The Stokes regions are 
determined by the Liouville transformation. 



Also we can deduce the last assertion of the theorem from |T^, §7.4]. This is proved in 113 



Theorem 5] for more general equations. □ 

Remark 1. Under the Liouville transformation, a neighborhood of infinity in each Stokes 
region in the complex 2;-plane maps to a neighborhood of infinity in either the upper or lower 
half Z-plane. So if u decays in a Stokes region Sj for some j, then u must blow up in the 
Stokes regions Sj+i and Sj-i. Otherwise, there would be a solution of the sine equation in 
the Z-plane which decays to zero in all directions. This is a contradiction. However, u might 
blow up in many consecutive Stokes regions (even in all Stokes regions) (see [|1^, §7.4]). 
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Definition. Let A G C and let u{z) ^ be an analytic function on C that satisfies (|l]). We 
say u is an eigenf unction and A is an eigenvalue, for if u{z) decays to zero along rays to 
infinity in the left- and right-hand Stokes regions (that is, if u has decaying asymptotics in 
(0), in these two regions). 

Remark 2. Given a Stokes region 5*^, there always exists a solution of —u"{z) + [Piz"^) — 
= Xu{z) that blows up in 5*^ §7.4]. So if there were two linearly indepen- 
dent eigenfunctions with the same eigenvalue, then all the solutions of —u"{z) + [P{z'^) — 
{izY^^^]u{z) = Xu{z) would satisfy m(±oo + Oi) = and there would be no solutions that 
blow up in the left- and right-hand Stokes regions. Thus there are no repeated eigenvalues, 
and all eigenvalues are simple. 



Remark 3. Note that if u{z) is an eigenfunction with eigenvalue A, then ui—z) is an eigen- 
function with eigenvalue A (an upper bar denotes the complex conjugate). So if an eigenvalue 
is real then u{z) = cu{—z) by Remark]^, and clearly |c| = 1. Writing c = e"^*"^ and replacing 
u by e'^'^u, we get that eigenfunctions with real eigenvalues are symmetric with respect to 
the imaginary axis. 

The main result of this paper is: 

Theorem 2. If X is an eigenvalue of then A 7^ and \ argA| < g;^- 

That the eigenvalues have positive real part was known already [|1^] (according to Mez- 
incescu |1^); our proof below includes a very simple argument for this fact. In the proof 
and elsewhere, we will use the following: 

Since u{z) decays exponentially along rays to infinity in the left- and right-hand Stokes 
regions, so does u' by the Cauchy integral formula. Therefore p(r) |u(re*^) p and p(r) |u'(re*^) p 
are integrable along the line r \—>- re*^ in C for any polynomial p{r), provided 16*1 < 2(2n+3) 
(so that the ends of the line stay in the left- and right-hand Stokes regions). 

Proof of Theorem ^. Let u be an eigenfunction with eigenvalue A, so that 

u"iz) + [-Piz^) + (2^)'"+X^) = -Xuiz), 

where P{z) = Til^^akZ^ for some afc > 0, A; = 0, 1, 2, n. 
Write 

A = a + z/5, a,/5e]R. 

Fix e with 1^1 < 2(2^- Let v{r) = u{re'^). Then v'{r) = v!{re^^)e'^ and v"{r) = u"{re'^)e^'^ . 
Thus our ODE becomes 

y"^r) + {[a + tf3- P{r\^'')y + ,2n+1^2n+lg^(2n+3)e| ^^^^ ^ ^ 
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Then we multiply this by e^**-^"'^^-'^w(r), integrate and use integration by parts to get 

(3) e-*(2n+3)e / |^,|2^^ 



/oo /"OO n 

'OO J —oo J — 



oo 



oo 



for all 16*1 < 2(2n+3) ' where we note that the line re*^ stays in the left- and right-hand Stokes 
regions where u (and hence u') decays exponentially to zero as z approaches infinity. 
Taking the real part of (||) gives (since |^| < 2(2n+3) ) 

/oo 
\v'\^dr 
■oo 

/oo 
\v\^dr 
■oo 

/oo 
Re [e-*(2n+l)ep^^2g2 
-oo 



If I dr. 



But Re [e-^(2"+i)ep(^2g2ie)j ^ T^^^akr'^^ cos{2n - 2A; + 1)^ > if > and \e\ < 
(certainly true if \6\ < 2(2n+3) )• '^'^ from (|^) we conclude that 

a cos(2n + 1)^ + /5 sin(2n + 1)^ > 0, 

a > \/3\tan{2n + l)e, 



for all \6\ < 2(2n+3) • That is 



for all < ^ < 

Taking 6* = gives a > 0, in particular A 7^ and taking 9 2{2n+3) S^^^^ 

a > /? tan^- '—. 

- 2(2n + 3) 

Then finally using tan0 = cot(^ — 4>), we have 

tan > — . 

2n + 3 a 

That is, |argA| < ^ 

Remark 4. We can extend Theorem |^ by allowing P to have some negative coefficients as 
long as P satisfies Re [e~*(2n.+i)eip(^^2g2i6'-jj ^ g |^| ^ Por example, with n = 3 and 

c G M, let P{z) = + cz^ + z; then Re [e-'^*^P(r2e2*^)] = [r"^ cos 6^ cr^ cos(36') + cos(56')]. 
So if cos2(3^) - 4 cos e cos(5^) < for |^| < f^, i.e. |c| < cos(f^) cos(ff ) ^ 1.837, then 



Re {e-^'^P{r'^e^^)\ > 0. So the theorem holds for this P provided c > -^Jf cos(f|) cos(ff ). 

Also by simple change of variables, we get the same result for H = — ^+[P(2;^) — (7(^2;)^"+^] 
for any non-zero real g. 

Moreover, by translations in C, we have the same result for H = + P{{z - 0^) - 
gj2n+i^^ — ^Y^'^^ for any ^ G C. For example, if u solves u"{z) — iz^u{z) = —Xu{z), then 
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v{z) = u{z + ai) solves v"{z) + [(3az^ — a^) —izlz"^ — 3a'^)]v{z) = —Xv{z) for any real number 
a. Observe v still satisfies the boundary conditions: f (±00 + Oi) = u{±oo + ai) = 0. 



Remark 5. The readers should notice that our boundary conditions are different, for n > 2, 
from those Bender and Boettcher ||^ take. In the zero boundary conditions of the 
problems —u" — {iz)^u = Xu for N > A are taken not on Stokes regions containing the real 
axis but instead on Stokes regions which are near the negative imaginary axis for large N. 

The next theorem extends Theorem 0. 

Theorem 3. Let A G C and n > 1. Suppose that u solves the ODE 

(5) u" - [P{z'^) + izQ{z'^)]u = -\u, m(±cx) + Oi) = 0, 

for some real polynomials P{z) = S^^Qa^^;^ and Q{z) = YTl^^h^z^ with all nonnegative 
and with 6„ G M — {0}. If for all k < n the coefficients a^, bk satisfy 

sm\2n - 2k)9 f ^^^^^^^ zfn=landk = 

6 77: , 77^ TT, T\Z^k < < 2afcafc+i if n > 1 and k = 0,n - 1 

cos 2n — 2k + Iw cos 2n — 2k — 1)6 -f ^ ^ ^1^?^ o 

^ I ; V ; ^ akttk+i if n> 1 and 1 < k < n — 2 

atO = ^T^r-TTT, then \ argA| < 



2(2n+3)' l"^6^N — 2n+3 " 

For n = 3, the coefficients of 6^ in (^) are approximately 3.41, 0.74 and 0.14 for = 0, 1, 2, 
respectively. 

Theorem |^ contains Theorem ^, just by taking = for k = 0,1, ...,n — 1 (in which case 
(D is trivially satisfied). 

Proof. The main idea of the proof is the same as that of the proof of Theorem Even if 
the equation (H) is little different from the equation for (|^), Stokes regions for @ are the 
same as for (0) (if 6„ has the same sign as (—1)""+^) or else are rotated by 180° (if 6„ has 
the opposite sign). See Section 7.4 in |jl6[ for details. But in either case, the lines r 1— > re*^ 
with \6\ < 2(2n+3) within the left- and right-hand Stokes regions, where we impose the 
zero boundary conditions. And this gives the integrabilities in the proof. 

Let v{r) = u{re^^). Then like we derived (^ in the proof of Theorem we have 

/oo 
\v\^dr 
-00 

/CO 
{cos[(2n + 3)e]\v'\'^ + S^^oOfc cos[(2n - 2k + l)e]i 
-00 



+ T^lzlbk sin[(2n - 2k)ey^+^\v\^] dr, for \e\ < 



TT 



2(2n + 3)' 

Since > for all k, we get a > by letting 6^ = in (^. (This is true for any Q with real 
coefficients and degQ = n.) 
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If we find conditions on and bk such that 

(8) ^l=Qak cos[(2n - 2k + l)ey''\v\^ + T^Hlbk sin[(2n - 2k)ey^+^\v\'^ > 0, 

for every r G M and l^'l < 2(2n+3) ' then we have from (|^) with 9 ± 2(2n+3) ' ^^^^ 

(2n + l)7r ^ (2n + l)7r 
acos \ \ ±^sin \ \ > 0. 
2(2n + 3) ^ 2(2n + 3) " 

So then | argA| < as desired, hke in the proof of Theorem |^. 

When > 3, we can rewrite the expression in @ as |f p times 

{flo cos(2?2 + 1)6' + rbo sin 2n6' + r^y cos(2n - 1)6'} 

(9) +Sfc=?{y cos(2n - 2A: + 1)0 + r6fcsin(2n - 2k)e + cos(2n -2k- l)e}r^'' 

+ { ^ cos 30 + rbn-i sin 20 + r2a„ cos } r2"-2 . 

Now (P) is nonnegative if each quadratic in r has non-positive discriminant: 

6osin^2ri0-2aoaicos(2n + l)0cos(2n-l)0 < 0, 
6fcSin^(2n-2A;)0-afcafc+icos(2n-2A; + l)0cos(2n-2A;- 1)0 < for 1 < k < n - 2, 

sin^ 20 — 2a„_ia„ cos 30 cos < 0, 

which is (H). The coefficients of bl in (||) are all increasing functions of < < 2(2^+3)' ^^'^ 
so it suffices that (H) hold at = 2(2^+3)- 

Now when n = 1, 2, it is easy to see similarly that the theorem holds. This completes the 
proof. □ 

Remark 6. In (§), the sign of r^'^+^lwpiir is difficult to determine because r can be 
negative as well as positive. 

The conditions in Theorem ^ are sufficient but not necessary, as is clear from the proof. 

We used = ^ + ^ to get (|) from (|^). If we use Ofc = 6kak + (1 — Sk)ak for some 
< 5/c < 1, we will get new sufficient conditions for the theorem. 

3. The zero-free region for u and u' 

The results in the previous section are based on the eigenfunction u decaying to zero as 
z approaches infinity in the left- and the right-hand Stokes regions. So consideration of the 
finite zeros of u may be useful for further results on our eigenproblem. 

For the next two sections, we will suppose H = —-£2 — (iz)^. See Figure || for the 
asymptotic behavior of the eigenfunction u. 

In this section, we provide a zero-free region for the eigenfunction u of 

(10) u" - iz^u = -Xu with u{±oo + Oi) = 
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/ \ 
/ \ 
/ \ 
/ \ 

Figure 2. In this Figure, rays are arg^; = i^, f , ff , ^ and A "+" 

indicates that the eigenfunction is blowing up while a "— " indicates that the 
eigenfunction is decaying to zero as z approaches infinity. 

and for its derivative m'. And we give some answers on how zeros of the eigenfunction should 
be arranged in C. 

It is obvious that u and m' do not share a common zero. Otherwise, by (p!0D, all the 
derivatives of u and u itself would vanish at the zero, and so m = 0. 
The following lemma is needed for our argument. Recall A = a + 

Lemma 4. Let 2; : [c, d] ^ C 5e a smooth curve with z'{t) 7^ fort G [c, d]. Ifu solves 
then writing z(t) = x{t) + iy{t), 

(11) ReK«)i:g 

/d I'd 
x'\u^{z{t))\'^dt + J [x'Re{iz%t) - X) - y'lm{tz%t) - X)] \u{z{t))\'^dt, 

and 

(12) Im(n'^)i:g 

pd I'd 

= - y'\u^{z{t))\'^dt+ [y'Re{iz%t)~X) + x'lm{iz%t)-X)]\u{z{t))\'^dt. 

J c J c 

Hille calls this lemma the Green's transform [p!6| , §11.3], and he uses it to get information 
on zero-free regions of solutions of linear second order equations (mainly with coefficient 
functions that are real on the real line). 

Proof. Let f{t) = u{z{t)) for t E [c,d]. Then f'{t) = z' {t)u' {z{t)) and 

= z'{t)u"{z{t)) = z'mz'it) - X]f{t). 
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A. 




Figure 3. Here Re {iz^ — A) is negative in since a > 0, while it is positive 
in Ai, A2 and A3. 



Hence by integration by parts, 



m 

z'it) 



fit) 



d I fi\2 



^dt+ I z'[iz^ - \]\f\''dt. 



Now by the formula f'{t) = z'{t)u'{z{t)) and splitting real and imaginary parts of the above, 
we get the lemma. □ 

Now we examine the consequences of the lemma. First, if Re {u'u) were not one-to-one 
on the imaginary axis, that would imply that the eigenvalue would be real by ( pT]) with 
z{t) = it. 

Remark 7. Second, another immediate consequence of Lemma ^ is that on any vertical line 
segments on which Im {iz^ — A) doesn't change its sign. Re ( function of y is one- 

to-one. On horizontal line segments on which Im {iz^ — A) doesn't change its sign, Im (Uxu) 
as a function of x is one-to-one (Mezincescu []TB|, §3] observed this last fact on the real axis, 
where y = 0). These observations are special cases of Hille's Theorem 11.3.3 in |l( 



Third, let us define open regions Aj and Bj, j = 1,2,3,4 as in Figure |^ and The 



following two theorems provide a large zero-free region for an eigenfunction u of ([T0|) and its 
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Figure 4. The level curves Im {iz^ - A) = with [3 > fixed. Here 
Im {iz^ — A) is negative in B4 while it is positive in Bi, B2 and B^. 




Figure 5. Im (-u'-u) < in the shaded area. Here di : Ivsiz = vW2, £2 
Im^ = -^/W/^ and £3 : Rez = -^/JJ^. 



derivative m', assuming A is non-real. Perhaps these theorems might help show that A must 
actually be real. The underlying ideas of the proofs are taken from Hille's book [ITB, §11.3]. 



Theorem 5. := ImA > then \m.{u'u) < on 



BiU {z e B^-.Rez < -\fWl} U {2 ^ : Imz > -^fjj2}, see Figure |. 
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Figure 6. Re (u'u) > in Figure 7. Re (u'u) < in 

the shaded area. the shaded area. 

Mezincescu §3] has previously observed that the eigenfunction has no zeros on the 
real axis, which obviously lies in the shaded region of Figure ^ Moreover, we see that all 
the zeros of u and u' in Imz > must be in Ai. 

Note that the lowest point in the closure d{B2) of B2 is —^(3/2 + i^(3/2. 

Proof of Theorem |^. For any ?/ G M, by (p!2[) with z{t) = t + iy and by ^(±00 + iy) = = 
u'{±oo + iy), we have that 

/X 
Im {iz{t)^ — A)!^!^*^^, 
-00 

and this is negative for x + iy G B4 with \y\ < ^/5/2, because then z(t) = t + iy G B4 for 
all —00 < t < X and so Im {iz(t)^ — A) < 0. 

This argument also shows that Im (u'u) < in {2; G i?4 : Rez < —■\/f3/2}; see Figure ^. 

Similarly in Bi, for all y we have that Im (u'u) = — Im (izit)'^ — X)\u(t + iy)\'^dt < 0. 

For z ^ Ai with Im^; > \f^j2 (so that z G A4), we use (|T^ along vertical line segments 
starting from points on the line Im^; = \f]5j2 to conclude Im (vlu) < in this region. □ 

Note that Im (u'm) = —\^\u(x ^ iy)'^ . So in the region in Theorem |, + iy)\ is an 
increasing function of y. 

Theorem 6. Assume (3 > 0. Then 

(i) Ke (u'u) > on the union of the regions A2, the region below A2 and the region TZ d B^ 
between A2 and B2 with the real part less than or equal to that of the zero uj^ of iz^ — X in 
the third quadrant. See Figure 

(a) Re (m'u) < on the union of the regions A^, the region below ^3, and the region in Bi 
with the real part greater than or equal to that of the zero 004 ofiz^ — X in the fourth quadrant. 
See Figure 0. 
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Obviously iz^ — A has three zeros. When (3 > 0, one of the zeros is in the second quadrant, 
one 0^3 in the third and one UJ4 in the fourth quadrant. Certainly these are the three points 
at which the boundaries of the Ai and Bi intersect. 

Theorem || with n = 1 shows that ^ < tan | ^ 0.73. It is easy to see that the rightmost 
point of 0/(^2) is at which lm{iz^-X) = x^-3xy'^-/3 > by < f3 < a. Thus 

the rightmost point of cl{A2) lies inside B3 as shown in Figure ^. Similarly, the leftmost 
point of 0/(^3) is — i), at which Im {iz^ — X) = — Sxy"^ — (3 < 0. Thus the leftmost 

point of 0/(^3) lies outside Bi as shown in Figure |^ 

Proof of Theorem |^. In the regions A2 and ^3, we use (0) with horizontal lines to infinity 
to get the statements in parts {i) and (ii) of this theorem. In the region TZ between A2 and 
B2 with the real part less than or equal to that of the zero of iz^ — A in the third quadrant 
(see Figure |^), we use ( |Tll) with vertical lines z{t) = x + it to show that Re (u'u) > 0. That 
is, we use Re{u'u)\^'^^^ = — lm{iz{t)'^ — X)\u{x + it)\'^dt. If x + id G TZ, we can find 
X + ic E d{A2 n -B3), so that Re [u'{x + ic)u{x + ic)] > 0, and — Im {iz(tY — A) > 0. Hence, 
the above integral is an increasing function of d. Hence we have the desired result in this 
region 71. 

The region below A2 is contained in B3 since the rightmost point of d{A2) lies in B3 (see 
Figure So a similar argument shows that Re {u'u) > in the region below A2. Also, the 
region below A3 is contained in B4 (see Figure 0) and so modified arguments show that the 
other statements of this theorem in part {ii) are true. □ 

Corollary 7. When ImX = P > 0, the zero-free region of u and u' contains the union of 
the three shaded regions in Figure and 

Note that in case Im A = /3 < we can get similar theorems corresponding to the above 
two, since u{—z) is an eigenfunction with eigenvalue A. The regions involved are simply the 
refiections of the above with respect to the imaginary axis. 

In case /5 = 0, so that A is real and A = a > 0, the regions Bi, B2, B3 degenerate to the 
sectors { — | < arg2; < {| < arg2; < ^}, { — < arg^ < — |} respectively, and we get 
the following theorem on zero-free regions. 

Theorem 8. Suppose X is real. Then lm{u'u) < on { — | < argz < ^} — Ai (which is 
a degenerate case of Figure while Re {u'u) behaves as in Figure and with Bi, B2, B3 
being sectors as above. 

Also Re {u'u) <0 m cl{Ai) n {Rez < 0} and Re {u'u) > in cl{Ai) n {Rez > 0}. 
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Figure 8. When A is real, the shaded area is the zero-free region of u and u'. 

Corollary 9. When A is real, the zero-free region of u and u' contains the union of all 
regions in Theorem see Figure |^. That is, u and u' can only have zeros in {iy : y > v^} 



U 



1^ G ^4 : -y < aigz < -^,lmz> -\/\/2^ U : |Re;z| < ^/\/2,lmz < -0^^ 



Remark 8. Bender et al. find numerically that u has some zeros along an "arch" within 
the unshaded region in Figure ^, when A is real. 

In proving Theorem we will use the following lemma. 

Lemma 10. Suppose C e Ai, Re (u'u) = at C- Then Re (u'u) < at all C - t E Ai,t > 
and Re (u'u) > at all ( + t ^ Ai,t > 0. 

Note that there is no restriction on the sign of ImA, in this lemma. 



Proof. On horizontal line segments z{t) = t + iy in Ai, ( [111) becomes 



Re(n'n)|;+:^= / \uMt))\^dt+ I Re {iz^{t) - X)\u{z{t))\^dt. 



I'd I'd 

I c+iy 

Since Re {iz^{t)—X) > in Ai, Re {u'u) is a strictly increasing function of x on each horizontal 
line segment in Ai. □ 

Proof of Theorem |^. The proofs of Theorems |^ and ^ give everything except the last state- 
ment of the theorem. For that, recall we can take u{z) = u{—z) by Remark ^ this implies 



14 K. C. SHIN 

Re (u'u) is an odd function with respect to reflection in the imaginary axis, so Re {u'u) = 



on the whole imaginary axis. Now we use Lemma 10 to complete the proof. □ 



By the last statement of Lemma |I], in the sector S^i^^^ that contains the negative imaginary 
axis, the eigenfunction u has only finitely many zeros. Now with the zero-free region in the 
Theorems ^ and ^, we see that u has only finitely many zeros in Im^; < 0. Since u has 
infinitely many zeros, u must have infinitely many zeros in Imz > 0. When (3 > (hence 
when /? < as well), by Theorem |^, u must have infinitely many zeros in Ai. Also when 
/? = 0, by Theorem ^ u has infinitely many zeros on the positive imaginary axis. 

The next theorem gives some information on how zeros of u and u' in Ai should be 
arranged, when f3 > 0. Note that all the zeros of u and m' in Im^; > lie in Ai by Theorem 
|. 



Theorem 11. Suppose u{z) is an eigenfunction of ( [7^ j with eigenvalue A G C with 1mA = 
/? > 0. Then 

(i) Ke {u'u) > for some point on the imaginary axis if and only if uu' has infinitely many 
zeros in Ai fl B2 and at most finitely many zeros in Ai fl -83 / and 

(a) Ke{u'u) < for every point on the imaginary axis if and only if uu' has no zeros in 
{z & Ai : lie z < 0} and infinitely many in {z & Ai : Re^; > 0}. 



We will use the following lemma along with Lemma pH 



Lemma 12. Assume ImA = /3 > 0. Suppose ReCi < ReC2 and Re (u'u) = at Ci, (2 
{where (i 7^ (2)- Then 

(i) G d{Ai n B2) =^ ReCi < ReC2 and lm(i < Im^2; and 
(a) C17C2 e cl{Ai n B^) => Red < Re(2 and lm(i > lm(2- 

Proof of part (i). We will first prove this for Ci; C2 G fl -82- Suppose that ReCi = ReC2- 
Then we could find a vertical line segment z{t) in Ai fl B2 whose end points are C,i and ^2- 
We apply (pT]) to this line segment to get 



= -y \m{iz'^{t)-X)\u{z{t))\^dt. 

This would imply -u = on the curve z{t) since Im {iz'^{t) — A) > in B2. So then since u is 
analytic, u = in C. This is a contradiction. Hence Re^i < ReC2- 

Similarly, suppose that Im^i > ImC2- Then we could find a smooth curve z{t) = x{t) + 
iy{t) in Ai fl B2 such that z{c) = Ci, z{d) = (2, x'{t) > and y'{t) < 0. Note that 
Im {iz^{t) - A) > and Re {iz^{t) - A) > in Ai fl i?2- This contradicts like for the case 
of ReCi = ReC2- We now see that the above argument still holds for Ci, C2 G d{Ai fl i?2)- 
Proof of part (ii). We use (^) again and a similar argument like in the proof of part {%). □ 
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Proof of Theorem |I7]. Suppose u{z) is an eigenfunction of (|TU]) with eigenvalue A G C with 
ImA = (3 > 0. Since u has infinitely many zeros in Ai (by the paragraph shortly before 



Theorem [Tl| ), certainly uu' also has infinitely many zeros in Ai. 

Proof of part (i). Suppose that Re (m'm) > for some point iy on the imaginary axis. By 
(pA|) with z(t) = it, we have that 

/•d 

(13) 



Re(u'w)|2 = /3 / \u{zy)\'dy 



So then since /? > 0, Re (u'u) > at every point it for t > y. Now by Lemma 10, we have 
that Re (u'u) > at every x + it E Ai for t > y and a; > 0. Thus uu' does not have any 
zeros in {z E Ai n B2 : Re z > 0,lm z > y}. 

The entire function uu' does not have infinitely many zeros in any bounded region. So 
if uu' had infinitely many zeros in Ai fl then uu' would have infinitely many zeros in 
{z E Ai n B2 : Re z < 0}. But if uu' has a zero zi in {z E Ai n B2 : Re z < 0}, then by 
(a), uu' has no zeros in {z E AiH B2 : Re zi < Re z < 0} . So then uu' would have 



Lemma 12 



infinitely many zeros in a bounded region. This is a contradiction. Thus uu' has infinitely 
many zeros in Ai fl B2 and at most finitely many zeros in Ai fl -Bg- 

Conversely, suppose that uu' has infinitely many zeros in Ai fl B2 and at most finitely 
many zeros in Ai fl i?2- Choose a zero zq in Ai fl B2. Then by Lemma we see that 
Re (u'u) > at ilm.zo since Re^o < 0. 

Proof of part (ii). Suppose that Re {u'u) < for every point on the imaginary axis. Then 
by Lemma 0, Re {u'u) < for every point in {z E Ai : Re z < 0}. So then uu' has no zeros 
in {z E Ai : Re z < 0}. Now since we know that uu' has infinitely many zeros in Ai, uu' 
must have infinitely many zeros in {z E Ai : Re z > 0}. 

Conversely, suppose Re {u'u) > for some point on the imaginary axis. Then uu' would 
have at most finitely many zeros in {z E Ai : Rez > 0} by the argument as in the proof of 
part {i). This completes the proof. □ 
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Remark 9. Since the negative imaginary axis is in the middle of a blowing-up Stokes region 
(see Figure |]), u{iy) blows up as y tends to — oo. On the other hand, the positive imaginary 
axis is a critical ray. We can show that \u{iy)\'^ < {const.)y~i for all y near positive infinity, 



by Theorem 7.4.4 in [p^. 



So the right-hand side of (|1^) approaches +oo as c tends to — oo (while d is fixed). Thus 
we see that Re [u'(ic)u(ic)] < for all c near negative infinity. However, the right-hand side 
of (O) is convergent as d tends to +oo (while c is fixed). So Re [u'u) may or may not become 



positive near infinity along the positive imaginary axis. 

The next lemma gives some information on zeros of u and u' in Im^; < 0, if any exist. 



There can only be finitely many such zeros, by the paragraph shortly before Theorem 111 



Lemma 13. Assume ImA = /? > 0. Suppose Im(^i < Im(^2 (ind Im (-u'-u) = Q at C^i, (2 
{where d 7^ ^2)- Then 

('>') C1X2 £ cl{A^ n B^) =^ IniCi < IniC2 and Re^i < ReC2; and 
(ii) C1X2 e cl{A4 n B4) =^ ImCi < ImC2 and ReCi > ReC2- 



Proof. We omit the proof because it is very similar to the proof of Lemma |T^. We use (|T2| 



instead of ([Tl|) , and also make use of Figures | and |^. □ 

Roughly speaking, then, the zeros move up and to the right in the third quadrant, and 
down and to the right in the fourth quadrant. This observation supports that when A is real, 
zeros of M in Im 2 < lie on an arch-shaped curve as in Figures 5 and 6 in ^ . 

4. Other properties of eigenfunctions 

Here we present a possible way of proving the conjecture that the eigenvalues X of H = 
— ^ — (izY are positive real. Given an eigenfunction u with eigenvalue A, Theorem |T^ 



below gives a class O of polynomials p{x,y) which are orthogonal to |up in the sense that 
f^^p{x,y)\u{x + iy)\'^dx = for all y. One can perhaps prove the conjecture as follows. 
Suppose Im A 7^ 0; if (9 is large enough then \u\'^ = 0, giving a contradiction. 
Let u be an eigenfunction oi H = —-^ — {izY with eigenvalue \ = a + i(3. 

Theorem 14. Let O = {polynomials p{- , ■) : f^^p{x,y)\u{x + iy)\'^dx = for ally}. Then: 
(i) x^ - 3xy^ -P eO, 
(a) for all m > 0, 

3y^^-/5^)(x3-3y^x-/5) 



m-\-l m + b m + 3 m + 2' 

-m(m - l)x'"-2 _ 4a;™(3a;2^ _ ^3 ^ ^) _ _il_^^m+2 ^ 

(in) if p G O then py + 2{x^ — 3xy'^ — (3) p{t,y)dt e O, and 

(iv) if p G O then p^x + Pyy + I2x'^yp + 4(x^ — 3x?/^ — /?) Py{t, y)dt G O. 
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For example the following polynomials are in O: 

P3ix,y) = - 3x?/^ - /?, by (i), 

pj{x,y) = x^-9a;V-5a;^/5 + 18a;y + 18a;2?/2/3 + 4x(/32-3?/), 

by applying (iii) to ps and multiplying by 2, 

Psix, y) = 2x^ - Wx^y^ - 7x^(3 + 30a;^t/^ + 2hx^y'^(3 + hx^{(3'^ - I2y) + lOy^ - 10a, 

5 

by applying {ii) with m = and multiplying by -, 
p^{x, y) = 2x^ - Ibx'^y^ - 6x^P + 27x^y^ + Ax^{P^ - 27y) + 2Axy^ + 2lx^y'^(3 - 2Axa, 

by applying {ii) with m = 1 and multiplying by 6, and 
Pio{x, y) = 20x^° - UAx^y^ - 55x^/3 + 252x%^ + l89x^y^P - 35x^(48y - (3'^) 

+420x^(y^ — a) — 210, by applying {ii) with m = 2 and multiplying by 105. 



We do not know whether Theorem 14 generates all the polynomials in O 



Proof of Theorem [I^. It is useful to have the following two formulas, which follow from 
multiplying (p!0D by u and separating real and imaginary parts: 

(14) Im [ux{x + iy)u{x + iy)]x = {x^ — 3xy'^ — P)\u{x + iy)]"^ 
and 

(15) Re [ux{x + iy)u{x + iy)]^ = + {-Sx^y + y^ - a)\u\'^. 

At the end of the proof we will justify the fact that we can differentiate through the integrals 
that follow. 

{i) This is clear by integrating (0), using the zero boundary conditions in the left- and 
right-hand Stokes regions. 

{ii) Suppose m is a nonnegative integer. Then 



7 poo l-OO O l-OO 

(16) — / x"'\u\^dx= x'^ — lul^dx = -2 x"'lm{uxu)dx 
dy i-oo J-oo oy 



2 



oo 

oo 



m + 1 



x''"~''^Im {uxu)xdx, 



where the last step is by integration by parts. 
So using dT^), we have that 



d f°° 2 f°° 

— x"'\u\^dx = / x"'+\x^ -3xy^ - (3)\u\^dx. 



dy J^^ m + 1 
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Hence, 



/ x'^\u(x + iy)\'^dx 
dy^ J-oo 



(17) 



m + 1 

Then again using the integration by parts and (0), we have that this equals 

2 



m + 



/CO ™m+5 „r?i+3 ™m+2 

[-6x-+^y + - 3y'^ - f3^){x' - 3xy' - (3)M'dx. 

171 + m + 3 m + 2 



Also, we differentiate ([T6|) without applying integration by parts: 



j2 /-oo 



(18) 



(19) 



/oo /»oo 
a;'"|upda; = -2 / a;'"[Re (m^.m)^ - 2|M^|2]da; 
-oo J —CXI 



-oo 

fOO 



/oo 
a;™[(-3x2?/ + - a)|u|2 - by (15). 

-oo 



Also, applying integration by parts twice to the right-hand side of ([18|), we have that ([T 
equals 



(20) 



/OO nc 
x'^''^\u\^dx + A I 
OO J —c 



X I 'Uq^ I diX • 



By equating (|^) and (pOf) , we get 



(21) 



m{m — 1) 



m-2i .|2 



X 



u\ dx 



X 



-Sx'^y + y^ — a)\u\'^dx. 



Hence equating ( P^D and (|I^) and substituting ( PT| ) give (ii). 
(m) Suppose that j'^^p{x,y)\u\'^dx = for all y. Then 



— / p{x,y)\u\^dx 
dy J-oo 



(22) 



[py{x,y)\u\^ - 2p(x,y)Im (M^M)](ix 



-oo 
oo 



Py{x,y)\u\ +2<j / p(t,t/)(it [>Im(M^M), 



by integration by parts. This with (0) gives {Hi), 
{iv) Suppose ^'^^p{x,y)\u\^dx = for all y. 
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Then we differentiate through (^) with respect to y again to get 

/oo 
[Pyy{.x,y)\u\^ - Apy{x,y)lm{uxu) - 2p{x,y)lvsi {iu^^u - iluxl"^)] dx 
-oo 

/oo r f px "s 

Pyy{x,y)\u\'^ + aIJ py{t,y)dt\lm{u:^u)x - 2p(a;, ?/)(Re (m^m)^. -2|m^.|^) 
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oo 
oo 



dx 



/ Pyit,y)dtUx - 3xy - /3) 



\u\ dx 



/oo 
pix,y) [-Re {uxu)x + 2{-3x'^y + y^ - a)\u\'^] dx, by (14) and 
-oo 



\u\ dx 



Pyy{x,y) + A^J py{t,y)dt\{x - 3xy - /3) - Ap{x,y){-3x y + y -a) 
Pa:{x,y)2Re {uxu)dx. 



But (y^ — a) /^p|-up(ix = 0, and so applying integration by parts again to the last term 
gives {iv). 

Now to complete the proof we need to show that we can differentiate through the above 
integrals, which reduces to showing that 



/ x"'\u\^dx 



dy 

So we estimate the following: 

r>00 -1 



d 
dy 



for each m > 0. 



d 



x'"^\-{\u{x + i{y + h))Y - \u{x + iy)Y} - 7— |u(a; + iy)\^]dx 



< 



—00 
00 



dy 
d 



r{l^(^ + ^(2/ + ^))l ~ Wi.^ + } - TrW{x + iy)\ \dx 



h 



< \h\ 



\x\ 



, dWu\ 

Qy2 



dy 



{x + i^{x))\dx, by the Mean Value Theorem, 



where — y\ < \h\. 

So then it is enough to show that for each M > 0, there exist Ci > and C2 > such 
that 



(23) \u{x + iy)\ + \u'{x + iy)\ + \u"{x + iy)\ < 



IC2I 



for all \y\ < M. 



exp[\x\ 

Proof of ( P5D is as follows: 

Suppose that z = re*^ with \9\ < ^. Then z is in the decaying Stokes regions (see Figure 
By the asymptotic expression (0) we get that for some C > 0, |-u(re*^)| < exp[— Irl*^] for 
all 1^1 < ^ and large |r|, say |r| > R. Also choose R > M + 1. 

Now the Cauchy integral formula says that 



u^^'Hz) = — 



"(C) 
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So then 

\u^'Hz)\ < — / HC)\\dC\ < fc!max{|w(C)| : |C - z| = 1} < A:!exp[-(|z| - 1)^] 

where the last inequahty holds if {C G C : |C--s| = 1} C {C G C : C = pe^<^, \p\ > R, |0| < 
Choose < C2 < C. Then exp[-{\z\ - 1)^] < exp[-\z\^^] if \z\ > R2 for large R2 > R. 
The region where 1^1 > R2 and {( e C : \C-z\ = 1} c {C G C : C = pe'^, \p\ > R, |0| < 

covers all of |?/| < M but a bounded region. 

Since the minimum of exp[— in this bounded region is strictly positive, we can find 

a large Ci > so that the left-hand side of (|23D is bounded by Ci exp[— l^;^^]. Thus ( p3|) 

holds since |x| < 1^1 . And this completes the proof. □ 

Corollary 15. Let u{z) be an eigenfunction of (IT^). Then \u{x + iy)\'^dx is a convex 
function. 

Proof. This is a consequence of (|20|) with m = 0, or it can be proved using the subharmonicity 
of ImP. □ 



5. Conclusions 

Using simple path integrations, we were able to prove that eigenvalues of (|l|) lie in the 
sector I argA| < 2^^^ and we extended the result for some more general Hamiltonians. Also 
we provide zero-free regions of eigenfunctions and their first derivatives, for the potential 
— {ix)^. Then finally we have the set O of polynomials p{x,y) which are orthogonal to |-up 
in the sense that f^^p{x,y)\u\'^dx = for all y. 

In a recent communication with Mezincescu, he pointed out that for the potential —{ix)^ if 
Im A = /3 7^ 0, combining < tan | with the equation (23) in [0 gives |A| > (| cos |)10^ ~ 
3 X 10^. So if any non-real eigenvalues exist, they are very large. 

In this paper we consider only polynomial potentials with odd degrees. However, a number 
of other authors have worked on even degree potentials, particularly quartic |TT], and 
sextic [|T], ^ polynomial potentials. Our techniques in proving Theorems ^ and |^ can be used 
to get information on eigenvalues for even degree potentials if both ends of a line passing 
through the origin stay in decaying Stokes regions. 

Obvious open problems are to narrow the eigenvalue sectors closer to the positive real 
axis, and finally to prove that the eigenvalues are real. Since some "PT-symmetric non- 
Hermitian Hamiltonians do not have all real eigenvalues, one might further want to classify 
PT-symmetric non-Hermitian Hamiltonians which do have positive real eigenvalues. 
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